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Abstract 

A  unified  formulation  of  radio  scattering  theory  in  the  language 
of  random  functions  is  given.  The  problem  of  probabilizing  atmospheric 
dielectric  noise  is  discussed  in  detail.  The  amplitude  of  the  received 
scattered  signal  is  a  complex  gaussian  random  function  of  time,  with 
mean  zsro.  Its  statistical  structure  is  therefore  completely  specified 
by  its  covariance,  which  is  simply  related  to  a  mixed  spectral  density 
of  the  scattering  dielectric  noise.  Local  properties  of  second-order 
random  functions  are  defined  and  discussed,  and  certain  gaps  in  the  liter- 
ature are  pointed  out. 
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1,  Introduction 

In  recent  year  it  has  been  discovered  that  the  scattering  of  radio  waves 
by  weak  fluctuations  of  dielectric  constant  in  the  earth's  atmosphere  (both 
troposphere  and  ionosphere)  can  serve  as  a  practical  means  of  radio  communi- 
cation. The  phrase  dielectric  noise  is  a  convenient  and  suggestive  designation 
for  these  fluctuations.  There  is  little  doubt  that  most  dielectric  noise  ih  the 
atmosphere  is  produced  by  the  action  of  turbulence  at  large  Reynolds  number, 
V.Tien  vre  wish  to  emphasize  its  turbulent  origin,  we  shall  refer  to  dielectric 
noise  as  dielectric  turbulence ,  borrowing  Eckart's  phrase'-  ^, 

The  present  paper  is  concerned  with  framing  the  radio  scattering  problem 
in  its  natural  language,  namely  the  theory  of  random  functions.  With  few  ex- 
ceptions, this  matter  has  received  only  a  desultory  treatment  in  the  literature. 
The  detailed  structure  of  atmospheric  dielectric  noise  is  still  very  much  in 
issue.  Although  for  the  purposes  of  this  paper  we  do  not  have  to  be  too  explicit 
about  the  microstructure  of  dielectric  noise,  we  assiune  that  it  is  produced  by 
turbulent  mixing  of  dielectric  constant  inhomogeneities,  and  is  described  by  the 
Kolmogoroff-Obukhoff  structure  law.  The  reasons  why  we  favor  this  model  of  di- 
electric turbulence  have  been  given  elsewhere  L  -l '  , 

2.  Preliminaries  on  random  functions 

The  natural  tool  for  probabilizing  the  radio  scattering  problem  is  the  theory 
of  random  functions,  treated  in  the  books  of  DoobL  -',  Loeve"-  ,  and  Bianc-Lapierre 
and  Forte t  L  -l,  to  which  the  reader  is  referred  for  many  details  and  occasional 
undefined  terms.  We  review  here  only  those  aspects  of  the  theory  which  are  in- 
dispensable for  our  present  purposes. 

The  locally  isotropic  pressure  field  intrinsic  to  turbulent  flow  is  a  negligible 
source  of  atmospheric  dielectric  noise,  as  shown  by  Villars  and  Weisskopf L  J , 
and  Gallet "-  ■' .  In  this  connection,  it  should  be  remarked  that  the  suggestion 
made  by  the  author'-  ,  that  dependence  between  Fourier  coefficients  of  the  velo- 
city field  may  greatly  enhance  turbulent  pressure  fluctuations,  is  fallacious. 
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In  what  follows,  r  denotes  the  position  vector,  referred  tc  an  arbitrary 
origin  fixed  with  respect  to  the  earth,  and  t  is  the  time.  If  for  every  value 
of  r  and  t,  the  quantity  f(r,t),  in  f^eneral  complex,  is  a  random  variable, 
f(r,t)  is  said  to  be  a  randon  function  or  ( random )  process,  of  space  and  time. 
Let  f(r,t,a)  be  an  arbitrary  realization  or  sample  function  of  the  process  f(r,t}, 
where  a  is  a  point  in  a  sample  space  S  which  indexes  the  realizations  of  f(r,t.;. 
The  statistical  structure  of  f(r,t)  is  specified  by  a  probability  measure  |j.,  de- 
fined on  the  sets  of  r,  or,  equivalently ,  by  the  joint  distribution  of  any  n 
random  variables  f(r  ,t  ),..,f(r  ,t  ).  Ensemble  averages,  so  called  to  dis- 
tinguish them  from  space  and  time  averages,  will  be  denoted  by  angular  parentheses. 
They  are  integrals  over  S  with  respect  to  the  measure  n,  or,  equivalently,  integrals 
with  respect  to  appropriate  distribution  functions,  derived  from  the  measure  p  in 
a  familiar  manner  (see  [6],  p.  10).  For  example 

<  f(r,t)  >  =  I  f(r,t,a)dti(c), 
S 


<  f(r,t)f*(r',t')  > 


f(r',t,a)f*(r'  ,t ',a)d^l(a), 


J 


where  the  asterisk  denotes  the  complex  conjugate. 

The  random  function  f(r,t)  is  said  to  be  second-order  if 

(1)  <  |f(r,t)|^  >  <  00  , 


for  every  r  and  t.     Equation   (1)   implies  that  the  mean     <   f(r,t)>  and   ttie   covariance 
<  f(r,t)f*(r  ',t  )  >  exist  for  every  r,  r  ,t,   and  t   ,     It  is  to  be  noted  that  in 
radio  scattering,   second-order  properties  of  the  received  signal,  e,g.,   its 
average  power,  its  envelope   correlation,  etc.,   are  of  primary  importance,  as  are 
the  second-order  properties  of  the  dielectric  noise,   specified  by  its  covariance. 


*  By  a  complex  random  variable   (function)  is  meant  a  pair  of  real   randon  variables 
(functions),  namely  its  real  and  imaginary  parts   (see    [7],  p.  152), 
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The  relevance  of  the  theorj'  of  second-order  random  functions  in  physical  problems 
like  radio  scattering  stems  from  a  natural  preoccupation  with  power-like  quantities. 
If  the  distribution  of  the  n  random  variables 

(2)  f(i\,   +  h,  t,   +  k),    ...,  fCr"    +  h ,  t     +  k) 

11'  'nn 

is  independent  of  n  smd  k,  for  all  r^ ,   t^ ,    ...,  r  ,   t   ,  h,  k,  and  n^   then  f(r,t) 
is  said  to  be   stationary  in  r  and  t,     Stationarity  in  r  and  t  S3p3rat3ly  is  defined 
in  the  obvious  way,  and  will  sometimes  be  called  r-ftationarity  and  t-stationarity. 
If  the  stipulation  made  in  connection  with   (2)   is  true   for  the  values  n  «  1  and 
n  ■  2,  without  reference  to  whether  it  is  true  for  all  n,   fCr'jt)   ii    said  to  be 
second-order  stationary  or  to  have  a  stationary  co variance »     Seco.id-order  station- 
arity is  the   case  of  most  interest  to  us  here,   in  keeping  with  our  preoccupation 
with  the  second-order  properties  of  f(r,t), 

Stationarity  as  just  defined  does  not  allow  for  the  trends  exhibited  by  most 
of  the  random  processes  that  occur  in  nature,     A  natural  way  to  take  such  trends 
into  account  is  to  introduce  the  weaker  concept  of  local  stationarity,  which  may 
be  qualitatively  defined  as  follows    .     Suppose  that  for  all  n,   r, ,   t,,,..,r  ,t   , 
the  distribution  of  the  n  random  variables   i^.)  is  independent  of  h  and  k  to  a 
good  approximation,  provided  that  the  points   (r, ,t, ),,.., (r  ,t   )  and  their  trans- 
lates  (r^+n,t^+k), ,.,,(r  +h,t  +k)  all  lie   in  a  certain  domain  D.     Then  D  is  said 
to  be  a  domain  of  local  stationarity  of  f(r,t),  and  f(r,t)   is  said  to  be  locally 
stationary  in  D.     Suppose  now  that  the  infinite  space-^ime  domain  can  be  divided 
into  non-overlapping  domains  of  local  stationarity  of  f(r,t),  and  let  D(r,t)  be 
the  domain  of  local  stationarity  containing  the  point  D(r,t),     This  state  of  affairs 
will  be   summarized  by  saying  that  f(r,t)  is  locally  stationary  in  the  (space-time) 

The  definition  is  qualitative  in  the   sense  that  it  contains  the  loose  phrase    "  to 
a  good  approximation"  ,  and  the  correspondingly  loose  concept  of    "  domain  of  local 

stationarity"  ,  but  is  adequate  for  the  purposes  of  this  paper.  It  is  patterned 

[91 
after  the  definition  of  local  homogeneity  given  by  Kolmogorof f "-  ■■ ,  in  his  first 

paper  on  the   statistical  theory  of  turbulence. 
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domains   n(r", t).     The  definitions  of  local  stationarity  in  r  and  t  are   apparent, 
and  involve,   correspondingly,    space  and  time  domains  of  local  stationarity.     In 
particular,  if  f(r,t)  is  locally  stationary  in  the  domains   D(r,t),  then  for  every 
fixed  t,  f(r,t)  is  locally  r-stationary  in  the  space  domains  T)At) ,   defined  as  the 
projections  on  the  x,y,z  axes  of  the  domains  D(r,t)j   similarly,   for  every  fixed  r, 
f(r,t)  is  locally  t-stationary  in  the  time  domains  D  (t),  defined  as  the  projections 
on  the  t  axis  of  the  domains  rt(r,t). 

Suppose  f(r,t)  is  locally  stationary  in  the  domains  D(r,t).     If  for  all   r  and 
t,  DiTft)  is  the  infinite   space-time  domain,  then  the  definition  of  local  station- 
arity reduces  to   the  ordinary  definition  of  stationarity,  which  we  may  call   plobal 
stationarity,   to  distinguish  it  from  local  stationarity.     For  n  =  1  and  n  =  2,  we 
have   a  definition  of  second-order  local  stationarity  or  of  locally  stationary  co- 
variance,   analogous  to  the  corresponding  definition  for  global  stationarity.     Again, 
it  is  the  case  of  second-order  local  stationarity  which  is  of  most   interest  to  us 
here.     Following  the  usual  heuristics,  we   shall  deal  with  a  random  function  f(r,t) 
which  has  a  locally  stationary  covariance  in  the  domains  D(r,t)  by  approximating 
it  by    "pieces"     of  suitable  random  functions  with  globally  stationary  co variances. 
The  piece  used  to  approximate  f(r,t)   in  the  domain  f(r,t)   is,  of  course,  taken 
from  a  globally  stationaiy  random  function  which  has  substantially  the  same  co- 
variance   in  D(r,t)  as  f(r,t).     The   qualitative  nature  of  this  approximation  scheme 
parallels  the  inherently  qualitative  nature  of  our  definition  of  local  stationarity, 
but  it  la  adequate  for  the  purposes  of  this  paper.     A  more  rigorous  treatment 
of  random  functions  with  locally  stationary  covariances  will  be  presented 
elsewhere. 
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From  a  given  random  function  f(r,t),  other  random  functions  can  be  derived  in 
a  variety  of  ways,  two  of  which  we  illustrate  here.  One  consists  in  assigning  a 
function  of  r  and  t  to  every  sample  function  f(r,t,a),  by  a  transformation 

^3)  T  f(r,t,a)  -  g(r,t,a)  ; 

we  then  define  a  random  function  of  space  and  time  g(r,t),  which  has  the  setmple 
functions  g(r,t,a),  and  the  same  sample  space  S  and  probability  measure  [i  as 
f(r,t).  The  other  way  consists  in  assigning  a  function  of  t  only  to  every  sample 
function  f(r,t,a)  by  a  transformation 

(U)  T  f(r,t,a)  -  g(t,a)s 

we  then  define  a  random  function  of  time  g(t),  which  has  the  sample  functions 
g(t,a)  and  the  same  sample  space  S  and  probability  measure  n  as  f(r,t). 

As  an  example  of  the  first  transformation,  we  center  the  sample  functions  of 
f (r,t)  about  their  mean.  As  an  example  of  the  second  transformation,  we  multiply 
each  sample  function  by  exp(iK»r)  and  integrate  the  product  with  respect  to  r  over 
a  finite  region  of  space.  As  we  shall  see  below,  this  example  is  a  special  case 
of  the  transformation  which  relates  the  dielectric  noise  responsible  for  radio 
scattering  to  the  received  scattered"  signal.  Of  special  importance  is  the  common 
case  where  the  transfonnation  T  is  homogeneous  with  respect  to  one  or  both  of  the 
variables  r  and  t,  either  everywhere  or  in  certain  domains  *  For  example,  if 
equation  (3)  implies  that 


T  f(r-r^,t-t^,a)  -  gir'T^,t-t^,a) 
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for  all  o,  r  ,   and  t  ,  then  T  is  said  to  be   (globally)  liomogeneous  with  respect 
to  r  and  t  (see    [B],  p.  19)*     If  equation  (U)  implies  that 

T  f(?,t-t^,a)     -     g(t-t^,a) 

for  all  a,   provided  t     lies  in  a  certain  domain  D( t)   containing  the  point  t,   then 
T  is  said  to  be  locally  homogeneous  with  respect  to  t  in  the  domain  D(t).     An  im- 
mediate consequence  of  the  homogeneity  property  is  that  one  can  infer  the  5ta^ion- 
arily   (global  or  local)  of  homogeneous  transforms  of  a  random  function  from  the 
slalionarity  of  the  random  function  itself,  as  illustrated  in  the  next  section. 
(For  a  stateroait  of  the  homogeneity  property  in  the  language  of  tiiue-invariant 
filters,  see  Laning  and  Battin,    [103 ,  p.  1^0.) 

In  technological  applications  of  the   theory  of  (globally)  stationary  raxidow 
functions,   it  is  necessary  at  some  point  to  invoke  an  ergodic  property   ,  i»e«,  to 
identify  ensemble  averages  with  suitable  experimentally  available  averages   (e.g., 
space  averages,   time  averages,   etc.).     Analogously,  a  local  ergodic  property  is 
needed  when  dealing  with  locally  stationary   random  functions,  but  one  must  be  ioore 
careful  about  aabuming  its  validity   (see  Sect.  2).     A  definition  of  such  a  property 
which  is  suitable  for  our  purposes  is  the  following.     Let  toe  random  function  f (x)  be 
second-order  locally  stationary  in  the  intervals  D(x),  where  T)(x)  denotes  the  domain  of 
local  stationarity  of  f(x)   containing  the  point  x.     (For  simplicity,  we  coniins 
ourselves  to  randoiu  functions  of  one  variable,  leaving  to  the   reader  the  straight- 
forward generalizations  for  random  functions  of  several  variables.)       Then  f(x)  Is 
locsLLly  ergodic  in  D(x)  if  the  following  two  statements  axe  true. 

1.     The  ensemble  average 
<  f(x)  > 
and  the  x-average 


; 


r 

f(x,a)dM.(a) 

S 


See,  for  example ,    (lO^*  P*  ■^3< 


-  7  - 


f(C,a)dC 

are  approximately  equal  for  approximately  all  sample  functions   ,     The   second 
integral  extends  over  the  domain  P(x),  of  length    X   [p(x)]].     Note  that  <  f(x)  > 
is  approximately  independent  of  x  in  D(x). 


2.       Let  x  and  x+x     be  two  points  in  the  same  domain  D(x)  =  D(x-t-x  )   such  that 
i.  QXx)]  »    |x    |.     Then  the  ensemble   average 

<  f(x)f*(x+x^)  >  =   J  f(x,a)f*(x+Xp,a)dti(a) 

S 

and  the  x-average 


f(£,a)f*(€+x^,a)dC 


are  approximately  equal  for  approximately  all  sample  functions.  Note  that 

<  f  (x; 

C(x). 


<  f(x)f  (x+x  )  >  is,  to  a  good  approximation,  a  function  of  x  only  in  the  domain 


Failure  to  define  and  use  ensemble  averages  has  led  to  some  minor  confusion 
in  the  radio  scattering  literature.  The  scattered  radiation  is  a  linear  transform 
of  the  dielectric  noise,  involving  an  integration  over  r  (see  equation  (16)).  Thus, 
if  the  correlation  function  of  the  dielectric  noise  is  defined  as  a  space  average, 
as  by  Booker  and  Gordon"-  -',  a  difficulty  arises  due  to  the  double  service  required 
of  the  variable  r.  This  difficulty  was  noted  by  Staras  •-  -1,  who  prefers  to  define 
the  correlation  function  as  a  time  average.  Since  no  integration  over  time  appears 
in  the  usual  approximation  for  the  scattered  signal,  this  may  be  considered  a  tech- 
nical improvement.  However,  it  is  easy  to  imagine  situations  where  time  integrals 

Jupt  as  almost  all  means  all  except  for  a  set  of  zero  n-measure,  approximately 
all  means  all  except  for  a  set  of  negligibly  small  ^-measure. 
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are  required,  e.g.,  if  the  bandwidth  of  the  incident  radiation  is  so  large  that  re- 
tardation effects  cannot  be  neglected.  One  way  to  avoid  these  bpurious  difficul- 
ties is  to  define  and  use  ensemble  averages,  and  to  invoke  ergodicity  only  at  the 
end  of  a  calculation,  when  it  is  finally  necessary  to  identify  ensemble  averages 
with  appropriate  space  and  tirne  averages.  On  the  other  hand,  it  cannot  be  denied 
that  there  are  versions  of  radio  scattering  theory  that  get  along  witnout  the  use 
of  ensemble  or  phase  averages. 

3«     Some  properties  of  dielectric  turbulence 

The   i^elative  dielectric  con&tant  in  the     atraotjphere  at  the  point  r  and  time  t 
is  a  second-order  random  function  6(r,t},  which  is  the  result  of  many  different 
kinds  of  dielectric  noise,  each  with  its  own  characteristic  periodicities.     The 
most  prevalent  kind  seems  to  be  dielectric  turbulence  at  large  Reynolds  nmober, 
which  in  this  paper  will  be  as  sums d  to  be  the  only  kind  of  dielectric  noise  present 
in  the  atmosphere.     Thus,  the  sample  functions  e(r,t,o),   the   sample  space  S,   and 
the  probability  measure  \i  will  be  taken  to  correspond  to  outcomes  of  the  process 
e(r,t)  that  differ  becau  -3  of  the  action  of  atmospheric  turbulence,   and  departures 
of  e(r,t)  from  etationaiity  will  be  ascribed  to  trends  (both  spatial  and  temporal) 
in  the  meterological  mechanisms  that  produce  the  turbulence. 

We  now  postulate  certain  properties  of  dielectric  turbiilence  that  are  sug- 
gested by  experiment.     We   admit  in  advance  that  the   complexities  of  meteorology 
ai'Q  such  that  the  best  that  can  be  hoped  for  is  that  these  postulates  are  an  ap- 
proximation to  the  true  state  of  affairs  most  of  the  time   and  in  large  regions  of 
the  atmosphere, 

1,  For  each  fixed  r,  c(r,t)  is  locally  t-stationary  in  time  domains  D  (t) 

of  the  order  of   several  minutes' .    (Airborne  refractometer  measurements  in  the 
troposphere  show  that  tliis  assmnption  is  usually  justified;  in  the  case  of   tne 
ionosphere  the  validity  of  the  assumption  is  more  hypotheticail,  aue  to  the  lack 
of  direct  measurements.)     On  the  other  hand,   6(r,t)  is  not  locally   r-stationary 
in  domains  large  enough  to  make  the  property  useful,  as  shown  by  the  spatial 
changes  in  <  e(r,t)  >  associated  with  the  normally  stratified  atmosphere,  not  to 
mention  the  abrupt  changes  associated  with  sharp  elevated  layers. 

2,  'ihe   cantered  second-order  random  function  A  e(r,t)  ■  6(r,t)  -  <   E(r,t)  > 

■  '  —    ^ 

is  locally  t-stationary  in  the   same  domains  D  (t)   as  e(r,t),  i.e.,  in  time  domains 

*As  we  shall  flee  below,  in  the  radio  scattering  problem  we  need  only   assume 
second-order  local  stationarity,  but  we  shall  not  always  state  this  explicitly. 
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of  the  order  of  several  minutes.  This  follows  from  the  fact  that  the  operation 
of  centering  is  locally  homogeneous  with  respect  to  t  in  the  domains  D  (t),  since 

<  e(r,t)  >  is  approximately  independent  of  t  in  D  (t).  Note  that  by  definition 

<  ^  e(r,t)  >  -  0. 

3.     For  each  fixed  t,     ^  e(r,t)  is  locally  r-stationary  in  space  domains 
D.  (r)  with  linear  dimensions  at  least  as  large  as  100  meters,    and  often  consider- 
ably larger.     The  validity  of  tl-ds  assumption  is  indicated  by  refractometer 
measurements  in  the  troposphere,  but,   again,   the   situation  in  the  ionosphere  is 
more  conjectural.     The  fact  that  ^  e(r,t)  is  locally  r-stationary  in  sizeable 
domains,  although  e(r,t)   is  not,   is  easily  explained  if  ^  e(r,t)  is  produced 
by  turbulent  mixing  operating  in  gradients  of  dielectric  constant.     In  this  case 
the  rms  dielectric  noise  power  <  ^  e   (r,t)  ->  ^     is  proportional  to  the  gradient 
of  the  mean  dielectric  constant  field"-  -I  *  L^-l   so  that 


(5)  <  A  e^(r,t)  >^/^    -  L^  grad  <   e(r,t)  >   . 

The   constant  of  proportionality  L     is  the  largest  eddy  size  of  the  turbulent  velo- 
city field     •     According  to  equation  (5),  the  presence  of  gradients  of  <  e(,r,t)  >, 

which  precludes  the  local  r-stationarity  of  6(r,t),   is  perfectly  compatible  with 
the  local  r-stationarity  of  ^  6(r,t),  provided  that  there   are   sizeable   domains 
in  which  grad  <  e(r,t)  >  is  substantially  constant,   and  in  which  the   velocity  field 
itself  is  locally  r-stationary,  as  seems  to  be  the  case.     The   space-time  domains 
in  which  ^  e(r,t)  is  locally  stationary    are  approximately  the  product   sets 

P(r,t)     -     D^(r)  X  Dj.(t)      . 


We  shall  use  the  symbol  "  to  mean  "  approximately  equal  to"  ,  sometimes  (as  here) 
only  to  within  an  order  of  magnitude;  -^   is  admittedly  a  symbol  with  a  somewhat 
variable  meaning,  but  this  disadvantage  is  offset  by  its  convenience  in  the  radio 
scattering  problem,  where  a  large  number  of  approximations  are  made. 


For  background  in  turbulence  theory,  see  G.K,  Batchelor 


&3] 


-lo- 
in radio  scattering  theory,  the  received  signal  is  attributed  to  the  centered 
random  function  /\^  e(r,t),  which  represents  turbulent  fluctuations  of  the  dielectric 
constant  about  its  mean   .     The  amplitude  of  the  current  in  the   receiving  antenna 
load  is  a  random  function  of  time  l(t),  which  is  a  weighted  Fourier  space -transform 
of  ^  e(r,t),   taken  over  a  finite   region  of   space  V  (see  equation   (16)),     Since 
the  transformation  l(t)   is  homogeneous  with  respect  to  time,  it  follows  that  l(t) 
is  locally  t-stationary  in  the  domains  D(t),  whose  dimensions    /(p(t)]   are  given 
by 


i[p(t)]       =     Min    J^  R(t)1    ,  r  in   V, 


where  D  (t)   are   the  domains  of  local  t-stationarity  of    A  e(r,t),  with  dimensions 
r 

j2.  [p  (t)].     If  all  the  domains  P  (t)  have  dimensions  of  the  order  of  a  few  minutes, 
then  so  do  the  domains  D(t),     The  local  t-stationarity  of  l(t)   in  domains  of  the 
order  of  a  few  minutes  is  an  experimental  fact.      (This  does  not  imply,  but  lends 
plausibility  to  the  local  t-stationarity  of  ^  e(r,t)  in  dom.ains  of  the  same   size.) 
For  example,   in  the  experiments  of  Chisholm  et  al  ^    -^ ,  the  envelope  distribution 
of  the  received  signal  was  found  to  be   approximately  Ray lei gh  if  four-minute 
samples  were  analyzed,  but  approximately  gaussian  if  16-minute   samples  were 
analyzed.     This  change  in  distribution  can  be  interpreted  as  a  slow  chanp^e  in  the 
parameter  of  the  Rayleigh  distribution,  which  describes  the  envelope  of  the  scattered 
signal  for  several  minutes  at  a  time. 


If  X  is  a  typical  point  of  a  set  E,  and  f  a  function  of  x,  then  the  minimua 
value  of  f  at  a  point  of  E,  and  the  maximum  value  of  f  at  a   point  of  E,  will  be 
denoted  by  Min  f(x),  x  in  E, and  Max  f(x),  x  in  E,  respectively. 
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li»     The  process   A  e(r,t)  is  characterized  by  a  correlation  distance  L  , 
which  is  the  linear  dimension  of  the  region  within  which  values  of  ^  e(r,t) 
are  appreciably  correlated.     As  the  symbol  L     indicates,  the   correlation  dis- 
tance is  a  way  of  characterizing  the  macro-eddy  dimension,  which  is  inherently 
not  defined  to  within  much  better  than  an  order  of  magnitude.     One  might,  for 
example,  estimate  L     as  the  distance  such  that  the  correlation  between  the  random 
variables  ^  6(r,t)  and  ^  e(r  ,t    )  is  <  ±flO,  provided  that    |r-r    |  >  L  •     Re- 
fractometer  measurements  in  the  troposphere  indicate  that  L     "^  100  meters  j 
estimates  of  L     in  the  ionosphere  vary,  but  the  value   L    ~  100  meters  is  favored 

O  r  ^    ,  Q 

by  some  authors,  e.g.,  Gallet L  ^»     In  analogy  with  the  situation  for  the  turbu- 
lent velocity  field  (see  below)  it  is  reasonable  to  afrume  that  the  joint  distri- 
bution of  the  random  variables  ^  e(r,t)  and  ^  e(r  ,t  )  is  approximately  gaussiar 
provided  that  |r-r  |  >  L  ,*  thus,  to  a  good  approximation,  ^  e(r,t)  and  A  e(r  ',t  ) 
are  independent  as  well  as  uncorrelated  if  |r-r  |  >  L^^*  By  a  correlation  volume 
is  meant  a  cube  of  side  L  , 

0 

5,     As   just  noted,   the  maximum  spatial  periodicity  (macro-eddy  size)  of 
A  e(r,t)  has  been  estimated  as  L    ~  100  meters,   in  both  the  ionosphere   aiid  tropo- 
sphere; the  corresponding  maximum  temporal  periodicity  (macro-eddy  lifetime)  is 

T     ~>    L  /v    'w    2  minutes,  where  v     is  the  rms  turbulent  velocity  (macro-eddy 
0  0      0  '  0  ^  -- 

velocity),   estimated  at  1  meter/sec'--'''-    J**.     Kolecular  dissipation  processes 
determine  corresponding  minimum  periodicities.     The  minimum  spatial  periodicity 
(  "  smallest  eddy"  )  of  dielectric  turbulence  is  approximately  the  same  as  that 
of  the  underlying  velocity  field.     It  has  been  estimated  as 

L  -^  Ji,  where    J?  is  the  parameter  in  the  correction  function  exp(- |r  1/ J? ),  used 
by  Booker  and  Gordon'-    -*»     Fstimates  of   H    are  summarized  by  Gordon  L    -' , 

Since  the  smallest  vinit  of  dielectric  turbulence  is  a  single  macro-eddy,   L     and 
T     are  estimated  lower  bounds  for  the  space  anu  time  dimensions  of  the  domains 


of  local  stationarity  of  ^  e(r,t). 
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Vj  ^  1  cm  (troposphere)'-  j*i-  'J 


^6) 

v|  ~  IC  meters  (ionosphere)'  , 

From  the  standpoint  of  the  radio  scattering  problem,  the  minimum  temporal  period- 
icity of  dielectric  turbulence  corresponds  to  the  case  where  the  mean  velocity 
field  <  V  >  convects  the  spatial  pattern  of  the  smallest  eddies  past  a  point 
fixed  with  respect  to  the  earth.  This  leads  to  the  estimate 

min    ' 

Since  the  mean  wind  velocity  certainly  never  exceeds  100  meters/sec  in  either 
the  ionosphere  or  the  troposphere,  we  find  with  the  use  of  (6)  that 

T  >  10   sec  (troposphere), 

T  .   i  lO'  sec  (ionosphere)  . 

min  ' 

If  <  V  >  <  V  ,  which  can  happen  only  if  the  turbulence  arises  from  causes  other 

than  shear  flow,  then  <  v  >  should  be  replaced  by  v  in  the  ex-ore ssion  for  T  .  , 

corresponding  to  the  convection  of  the  smallest  eddies  by  the  largest  eddies. 

This  results  in  estimates  of  T  .  which  are  about  IOC  times  lareer  than  the 

mm  "^ 

above  values.     These  considerations  are  enouph  to  show  that  the  fading  of  the 
scattered  signal  involves  both  the  lifetimes  of  the   scattering  eddies  and  their 
convection  by  the  mean  wind   or  by  the  macro-eddiec      . 


This  value  of    y}   is  obtained  bv  using  the  estimates  L     ~  100  meters  and  v    ~  1 
meter/sec;  y)  is   only  ''•reakly  dependent   on  L  ,   so  that  Villars  and  Weisskopf '- -^ 
arrive   at   the   same  estimate  of  Vi  ,   using  a  much  larger  value  of  L   .    A  value  of  vi 
of   this  order  of  magnitude   follows  from  an  analysis   of  the   scaled  radio   scattering 
experiments  of  Bailey,   Bateman,   and  Kirby '-    -' ,   as  pointed  out  by  ^lOieelon'-    -' » 


-  13  - 


6.     We  must  now  postulate  a  plausible  ergodic  property  in  order  to  identify 
ensemble  averages  of  A  6(r,t)  and  certain  of  its  transforms  with  suitable  space 
and  time   averages.     A  cautionary  remark  is  needed  before  invoking  a  local  ergodic 
property  like  that  described  in  Section  2,    since,   as  already  noted,  it  is  con- 
ceivable that  the  dimensions  of  some  of  the  domains  of  local  stationarity  of  the 
procesp  /\6(r,t)  may  be  of  the  same  order  of  magnitude  as  the  macro-eddy  siie  L^ 

and  the  wacro-eddy  lifetime   T  .     It  is  clear  that  sample  fxinctions  measured  in  such 

o 

a  domain,  with  dimensions  L  and  T  ,  will  not  be  large  enough  to  exhibit  typical 

'  0  0 

behavior  of  the  process  insofar  as  events  on  the  space-time  scale  of  the  macro- 
eddies  are  concerned.     To  avoid  this  difficulty  we  uaist  take  care  to  apply  the 
local  ergodic  property  only  to  transforms  of  ^e(r,t)  which  are  substantially  un- 
affected by  macro-eddy  behavior.     This  point  of  view  is  especially  cogent,  when  it 
is  recalled  that  description  of  events  on  the  space-time  scale  of  the  macro-eddies 
lies  outside  the  province  of  statistical  tuibulence  theory,  which  by  its  very  nature 
can  make  statements  only  about  the  local  structure  of  turbulence.      Kolmogoroff 
made  this  fact  explicit  in  his  original  researches  on  turbulence"-  -•   by  introducing 
the  structure  function 

D^{r')     -    <  |v(r,t)  -  v(?*  r',t)t^, 
where  v  denotes  some  component  of  the  turbulent  velocity  field.     Statements  about 
D  {?' )  are  then  made  only  for    jr '  |  «  L  ,  so  that  ths  effects  of  the  macro-eddies 
are  automatically  subtracted  out.     The   amplitude  I(t)  of  the  received  scattered 
current  is  another  important  example  of  a  trarisform  of  tiie  process  ^  e(r,t)  which 
is  insensitive  to  macro-eddy  events,  since  it  is  produced  by  scattering  eddies  at 
least  an  order  of  magnitude  smaller  than  L^.     Correspondingly,   the  largest  periodi- 
cities of  I(t)  are  much  less     than  T    «     With  these  qualifications,  we   shall  assume  a 
local  ergodic  property  of  the  type  discussed  in  Section  2,  whenever  needed. 

*In  fact,  it  is  often  difficult  to  distinguish  between  macro-eddy  events  and  the 
meteorological  events  responsible  for  the  turbulence. 
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7,     It  is  sometimes  assumed  that  ^  e(r,t)   is  a  gaussian  process.     There  is 
very  little  likelihood  that  this  is  the  case.     Since  the  turbulent  velocity  field 
is  not  a  gaussian  process   ,  other  turbulent  fields  should  be  non-gaussian  to  the 
extent  that  they  share  the   structure  of  the  underlying  velocity  field.     Therefore, 
it  is  reasonable  to  expect  that  atmospheric  dielectric  noise,  which  is  due  primari- 
ly  to   the  turbulent  mixing  of  dynamically  passive       scalar  quantities  will  be  found 
to  be  non-gaussian'      .     Presumably,   departures    from  the   gaussian  distribution 
will  first  show  up  in  the  bivariate  distribution  of  ^  e(r,t)   at  points  separated 
by  less  than  L  ,   in  analogy  with  the   situation   for  the  turbulent  velocity  field. 

0 

It  is  recommended  that  refractometer  data  be  analyzed  with  a  view  to  studying  the 
probability  distribution  of  ^e(r,t)',  this  does  not  seem  to  have  been  done  as 
yet. 

The  reader  will  note  that  our  model  represents  turbulent  fluctuations  about 
the  actual  mean  dielectric  constant  field  <  e(r,t)  >,  rather  than  about  the  vacuum 
value  of  unity,  as  often  written.  We  have  chosen  this  model  for  two  reasons. 


*  Although  the  turbulent  velocity  field  is  univariate  gaussian,  its  bivariate  dis- 
tribution is  markedly  non-gaussian  for  small  spatial  separations,  but  becomes 
gaussian  when  the  points  are  far  enough  apart  so  that  the  corresponding  two-point 
velocity  correlation  is  quite  small.  See  Batchelor"-  -' ,  chapter  0.  A  gaussian 
velocity  field  could  not  represent  turbulence,  since  the  vanishing  of  its  third 
moments  would  preclude  the  transfer  of  e nergy  from  larger  to  smaller  eddies  char- 
acteristic of  turbulent  flow.  The  non-gaussian  character  of  the  velocity  field 
corresponds  to  statistical  dependence  between  the  Fourier  coefficients  of  the  velo- 
city field, 
^*A  turbulent  field  quantity  like  temperature  is  said  to  be  dynamically  pastive  (in 
a  given  turbulent  regime)  if  any  turbulence  produced  by  its  inhomogeneities  is 
negligible  compared  to  the  turbulence  produced  by  other  causes. 
""As  shown  below,  departures  of  ^  e(r,t)  from  the  gaussian  distribution  do  not 
affect  the  scattered  signal.  Thus,  if  we  are  not  concerned  with  the  detailed 
structure  of  the  dielectric  turbulence,  we  can  replace  ^  e(r,t)  by  the  zero- 
mean  gaussian  process  with  the  same  covariance,  without  in  any  way  changing  the 
scattered  signal. 


A«-*M-  . 
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First  of  all,   since  we  favor  turbulent  mixing  operating  in  pre-existing  inhorno- 
geneities  of  dielectric  constant  as  the  source  of  most  atmospheric  dielectric 
noise,  it  is  desirable   tc  exhibit  the  terra  <  e(r,t)  >,  which  explicitly  represents 
the  inhomogeneities  of  dielectric  constant.     Secondly,  although  we  do  not  attri- 
bute the  received  signal  to  <  e(r,t)  >,   it  is  not  to  be  denied  that  the  mean 
dielectric  constant  field  has  an  effect  on  radio  scattering,   an  effect  convention- 
ally concealed  in  the  antenna  pattern  propagation  factors.     Writing  out  the  term 
<  e(r,t)  >  serves  as  a  reminder  that  one  must   justify  neglecting  the  departure 
of  <  e(r,t)  >  from  unity   at   any  stage  of  a  calculation,  especially  since  the 
departure  of  <  6(r,t)  >  from  unity  is  usually  much  greater  than  the  largest  fluc- 
tuations of  dielectric  constant  produced  by  turbulence, 

k.       Harmonic  analysis  of  ^  6(r,t) 

Of  particular  interest  in  physical  problems  is  the  representation  of  random 
functions  as  Fourier-Stieltjes  integrals.     Accordingly,  we  shall  assume  that 
/\  e(r,t)  is  harmonizable,   i.e.,  that  there  exists  a  second-order  random  function 
g(k,co),  with  a  covariance  of  bounded  variation,  such  that 


(7) 


^  e(r,t)     «         dg(k,w)  exp  -i(k   •  r  +  wt)      ,      (q.m.) 


where  all  variables  of  integration  range  from  -cd  to   oo,   as  they  will  in  any 

integral  that  we  write  without  limits.     Equation  (7)  is  a  stochastic  integral,   in 

the  sense  that  it  corresponds  to  the  limit  of  a  sequence  of  random  variables.     Let 

S  (r,t)  be  a  suitable  approximation  to  the   stochastic  integral  which  is  the  sum  of 
n 

For  a  discussion  of  hennonizability,  see  Loeve '- ^  p.  U7U,     Note  that  we  use  a 
minus  sign  in  the  complex  exponential  instead  of  the  more  usual  plus  sign.  Sto- 
chastic integrals  of  the  type  used  here,  i.e.,  Rieraann-Stieltjes  integrals  in 
the  quadratic  mean,   are  discussed  in   [7],  p,  U72. 
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n  random  variables.     The  abbreviation   (q.m.)   appearinp  in   (7)   signifies  that   the 
convergence  of  the   sequence   S   (r,t)    to  ^  e(r,t)   is  meant  in  the   sense  of  the 
quadratic  mean,   i.e.,  that 


lijt)       <    I  A  ^^r,t)   -  S^(r,t)|^>      -     0, 


for  every  r  and  t.     Of  the  reasons  for  this  stipulation,  perhaps  the  most  cogent 
is  that  the  integral   (7)  may  not  exist  for  the  individual  realizations  of 
g(k,co)".     It  follows  from  the  fact  that  ^  e(r,t)  is  a  real-valued  random 
function  that 

dg*(I^,co)     =     dg(-lf,-aj)j 

moreover,  <  dg(K,co)  >  =  0,  sance<  ^  e(r,t)  >  ■  0. 

Since  we  intend  to  approximate  A  e(r,t)  in  its  domains  of  local  stationarity 
by  pieces  of  suitable  globally  stationary  random  functions,  as  discussed  in  section 
2,  it  is  enough  to  consider  the  case  where  the  covariance  of  ^  t(r,t)  is  globally 
stationary  in  r  and  t.  Then,  the  covariance  of  ^  e(r,t)  is  a  function  only  of  the 
coordinate  differences  p  ■  r-r  and  L  =  t-t  ,  i.e., 

<  /\  e(r,t)  A  e(r',t')  >  ~  <  /\  t^  >   R(p,t), 

where  R(p,  X)  is  the  space-time  correlation  function  of  ^  e(r,t).     F^irthermore, 
by  a  well-knovm  theorem  of  probability  theory   (see    [7],  p.  U83),   the  covariance 
of  dg(?,'xO   is  of  the  form 

(8)  <  dg(Ic,co)dg*(k  ',co')  >  =  <  A  e^>  Y(k,a3)5(k-k')6(a)-m')dkdk"'dwdco'. 

where  the  non-negative   function  Y(k,cci)  is  the  space-time  spectral  density  of 

A  e(r",t),  measured  in  units  of  the  variance  <  ^  t     >,   so  that 

___  -  ^ 

If  the  process    A  e(r,t)   is  (plobally)  stationary,   almost  all  realizations  of  g(k,co) 

are  of  unbounded  variation,  so  that  almost  all  the  Riemann-Stieltjes  integrals  of 

the  individual  realizations  of  g(k,a))  fail  to  exist. 


-  17  - 


T(k,co)  dk  do)  »  1. 

(We  can  assume  that  the  process  ^  e(r,t)  has  no  line  spectrum,  an  assumption 
whicn  can  be  readily  justified  for  the  turbulent  velocity  field  (see  [l3J, 
pp»  25,  85),  and,  presumably,  for  dielectric  turbulence  as  irell.)  It  is 
easily  seen  that  R(p,'?')  and  Y(K,a))  are  Fourier- transform  pairs,  i.e.,  that 

R(P,  T)  =   T(k,to)  expUi(k  •  p  +  co  ir)|  die"  dco  , 

T(k,ai)  =  l/(2n)'^   R(p,r)  exp|i(k  •  p  +  w'C)1  dp^df. 

As  already  noted,  finite  Fourier  space-transforms  cf  the  process  ^  e(r,t) 
are  of  great  importance  in  radio  scattering  theory.  Accordingly,  we  introdace 
the  harmonic  decomDOsition 


(9) 


^  e(r,t)  -   dg^(k,t)  exp(-ik  •  r),    (q.m.) 


corresponding  to  spatial  harmonic  analysis  of  ^  e(r,t)  at  the  time  t.  The 
covariance  of  dg-(K,t)  is 

(10)        <  dg^(k,t)dg*(k ',t')  >  =  <  ^  e^  >  J  (k,t-t ')5(k-k ')dk  die  '  , 
where,  as  is  easily  verified 

^(k,X)  »  j  T(if,co)exp(-iojr  )dt  -  l/(2n)^  j  R(^,  r )exp(ik«p)dp 

The  mixed  quantity  ^  (k,  t)  is  a  spectral  density  in  the  variable  k,  and  an 
unnormaiJ.zed  correlation  function  in  the  variable  f  *  t-t  •  Consider  now 
the  finite  Fourier  space- transform 
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(11)  J(d,K;t)   «   J     /\  e(r,t)  exr(iK«r)   dr   , 

where  V  is  a  cube   of   side   d,   centered   at  the  origin;    J(d,K;t)   is  a  transform 
of    A  e(r,t)   of  the   type  of  equation  (U).      Substituting  (9)   in{ll),  we   obtain 

(12)  J{d,t',t)     =     V       sinc|(K-k)d/2J  dg^(lf,t),      (q.m.) 

where  sine  k  is  tho  function 

sine  k  =  (sin  k  /k  )(sin  k  /k  )(sin  k  /k  ) , 

Using  (10),  we  find  the  followinp  expression  for  the  covariance  of  j(d,K}t) 

(13)  <  J(d,K;t)J''"(d,K;t')  >  -  V^  <  ^  e^  >   sinc^  (K-k)d/2  ^  (k,1:)dk  , 


where  T  =  t-t  . 


The  function  sine  Ic  has  the  property  tbnt 


(Ih)  1  sine  Ic'  dlf 


2  -».  -X   3 

sine  k  dk  =  n 


Moreover,  most  of  the  contribution  to  the  integrals  (lU)  comes  from  a  region 

of  volume  n  in  k-space,  centered  about  the  origin.  Correspondingly,  most  of 

the  contribution  to  (13)  comes  from  the  volume  V  ■»  (2n)  /V  in  k-space,  centered 
about  the  point  if,  so  that  we  can  write 

(15)  <  J(d,K;t)/(d,Kit')>-  (2ti)^  V  <  /N^  e^  >  I  (K,r) 

to  a  pood  approximation,  provided  that  the  following  two  qualitative  conditions 
are  met.  (Note  that  as  d  ->oo,  the  approximation  (l5)  becomes  an  exact  equality.) 

1,  For  every  X   ,  ^(k,'^)  must  be  substantially  constant  in  the  volume  V 
centered  about  t.     For  any  reasonable  shape  of  ^(k,  tT),  this  means  that  K  +  2tT/d  ~  K, 
where  K  «  |  {^  | .  If  vre  divide  both  sides  of  this  approximate  equality  by  K,  it 
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becomes  2n/Kd  «  1,   or  Kd  »   2n.  This  condition,   i.e.,   that   the  linear  dimension 
of  the   finite  sample  of  ^  e(r,t)  be  long  compared  to  the  wavelength  at  which  the 
spectrum  is  being  estimated,   was  noted  by  Batchelor '-    -^ » 

2.     J{k,  t )  must  not  become  so  large  away  from  the  point  K  that  most  of  the 
contribution  to  the  covariance  of  J(d,lfjt)   comes  from  a  region  of  k-space  lying 
outside  the  volume  V     centered  about  K.     If  this  condition  is  insured  by  the  first 
condition,  namely  Kd  »  2n,  then  we   shall  say  that  ^  e(r,t)  has  the  local  harmonic 
property   (in  r),  or  that  its  spectral  density  has  the  local  harmonic  property  (in 
k).     It  should  be  noted  that  the  approximation   (15)  is  always  good  if  d  »  L  j   it 
is  assumed,  of  course,  that  there  are  values  of  d  «  L  ,  which  satisfy  the  rela- 
tion Kd  »  2iio    Since  ^kjT)  is  a  correlation  function  in  the  variable    t   ,   it 
follows  that  ^(k,  T)  <  ^1<,0),   so  that  if  the  space  spectral  density  ^(k,0)  has 
the  local  harmonic  property,   so   should  the    "  riixed*  spectral  density  5(k,  f). 
There  are,  of  course,  random  functions  without  the   local  harmonic  property   .     How- 
ever, it  is  intuitively  evident  that  turbulent  fields  in  general,   and  dielectric 
turbulence  in  particulai-,  have  the  local  harmonic  property.     For  in  the  picturesque 
but  meaningful  language  of  statistical  turbulence  theory,  harmonic  content  at  wave 
number  K  corresponds  to  the  presence  of  eddies  of  size    ~    2n/K,  which  are  only 
weakly  coupled  to  eddies  of  appreciably  different  sizes   ".     It  follows  that  a 
"piece"     of  turbulence  large  snoueh  to  contain  many  eddies  of  size  ~  2ti/K  should 
be  large  enough  to  determine  the   spectral  density  of  the  turbulence  about  wave 
number  K,  i.e.,  that  the  turbulence  has  the  local  hannonic  property.     Thus,  we 
shall  assume  that   (l5)  is  a  good  approximation  if     ^    :j  10,   say;  Batchelor  L    -' 
makes  the    same  assumption,  although  less  explicitly. 

For  example,  the  output  of  an  RC  filter  with  a  white  noise  input  does  not  have  the 

local  harmonic  property, 

A  basic  assumption  of  statistical  turbulence  theory  is  that  eddies  of  a  given  size 

"  forget"     the  eddies  of  much  larger  size  which  give  rise  to  themj  if  this  were  not 
the   case,   one  would  hardly  be    justified  in  constructing  a  statistical  or  proba- 
bilistic model   of  turbulence. 
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It  is  not  to  be  denied   that   these   heuristics  should  be  backed  up  by  some 
rigforous  mathematics.     The   author  has  begun  this  study,   and  plans  to  present  the 
results  in  a  future  paper. 

5.       Relation  between  the  scattered  signal  and  ^  e(r,t) 

Itie  geometrical  configuration  of  a  radio  scattering  experiment  consists  of  a 

transmitting  antenna  A  emitting  electromagnetic  radiation  that  is  scattered  into 

a  receiving  antenna  A     by  dielectric  noise  in  the  common  volume  of  the   antenna 

beams.     Before  writing  the  approximation  for  the  received  scattered  current,   ve 

Edve  our  versions  of  some  rather  standard  nomenclature  and  symbols. 

1,  c  is  the  velocity  of  light;  co  ,   the   angular  carrier  frequency;  k  »  u  /c, 
the  carrier  wave  number;    K  =  2nc/oo  ,  the  carrier  wavelength. 

2.  D  is  the  distance  from  A  to  A     along  the  earth's  surface   ;  R(r),   the  dis- 


•  /-A, 


tance  from  4  to  the  variable  point  rj  R  (r),  the  distance  from  A  to  r. 

3.  Q{T),gi{r)   are  the  polar  angles  of  r  referred  to  A  as  origin;  9  (r), 

^   (r),  the  polar  angles  of  r  referred  to  A  as  origin;  n(r),  the  unit  vector 


with  polar  angles  ©(r),5((r)  pointing  away  from  A;  n  (r),  the  unit  vector  with 
polar  angles  ©  (r),  0  (r)  pointing  toward  A  . 

U.  if(r)  is  the  scattering  wave  vector  associated  with  the  point  r,  i.e., 


the  quantity  ,^,r^.,^^   ^.^"T 

K(r)     =     (2tt/X)|  n  (r)  -  n(F)J    ; 

(5>(r)  is  the   scattering  angle  associated  with  t  he  point  r,   i.e.,  the  angle  between 
n(?)   and  n  (r).     For  the  usual  small  sc^ittering  angles,   the  magnitude  of  K(r)  can 
be  written  as 

K(r)   -    2ti  (a>(r)A  . 


We  assume  that  A'   is  well  beyond  the   geometrical  horizon  from  A,  i.e.,  that 
there  is  no  optical  path  between  A  and  A'.     The  effect  of  dielectric  noise  on 
line-of-sight  propagation  is  a  problem  requiring  its  own  special  methods,    and 
will  net  concern  us  here. 
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?      ?      I  I    „ 

5,  f(e,^)  and  f   (©,(?(   )   are  the  pattern  propagation  factors  of  A  and  A       j 

F(r),   the  joint  pattern  propagation  factor  associated  with  the  point  r,   defined 
by 

F(i^     =     f|©(r),  0(?)1       f '[©'(?),  g\t)]    . 

t 
We  assume  that  A  and  A  are  sited  so  that  their  directions  of  maximum  gain  inter- 

sect  in  a  point  r  ,  at  which  F(r  )  -  1,  We  also  assume  that  both  A  and  A  pro- 
c  c 

duce  linearly  polarized  radiation,  and  are  sited  with  their  polarization  directions 
parallel, 

6,  The  region  of  space  where  F(r)  is  substantially  different  from  zero  is 
called  the  (nominal)  scattering  volume  V,  and  the  point  r  is  called  the  center  of 
V.  Although  V  serves  to  indicate  the  region  of  space  responsible  for  the  received 
scattered  radiation,  most  of  the  radiation  may  come  from  a  subset  of  V,  known  as 
the  effective  scattering  volume  (see  below), 

7,  <S)  j^  is  the  minimum  scattering  angle.  Using  the  •  four- thirds  earth" 
approximation  to  estimate  the  effect  of  atmospheric  refraction  (see  [20],  p»  53), 
we  have  (5)  ^   "^  3r)/Ua,  where  a  is  the  radius  of  the  earth. 

We  now  write  an  approximation  for  the  signal  scattered  by  weak  dielectric 
noise  which  has  become  standard  since  its  introduction  by  Booker  and  Gordon^  -J , 
It  states  that  the  complex  amplitude  of  the  current  in  the  receiving  antenna  load, 
measured  in  units  of  the  amplitude  of  the  load  current  that  would  be  found  if  the 
antennas  were  in  free  space  and  oriented  with  their  directions  of  maximum  gain 
along  the  line  joining  them,  is  given,  to  a  good  approximation,  by 

(16)  I(t)  ~  (nD/X^)    --^, A'^(r,t)exp[iir(r)-r]  dr  } 

{  R(r)R  (?)  L     J 

By  the  pattern  propagation  factor  is  meant  the  quantity  defined  by  Kerr  (  [20], 
p.  28),  It  is  the  analog  of  the  antenna  pattern  function  ([20],p.3U),  and,  in 
principle,  incorporates  all  the  propagation  effects  of  the  actual  surroundings 
of  the  antenna,  including  reflections  from  the  earth's  surface,  and  refraction  in 
the  earth's  atmosphere.  In  practice,  rough  approximations  are  used,  due  to  the 
difficulty  of  ascertaining  precisely  the    'surroundings'   of  the   antenna. 
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the  integral  is  over  the  scattering  volvime  V,  and  an  irrelevant  phase  fsctor  has 
been  omitted.  The  interpretation  of  (16)  is  that  l(t)  is  a  complex  random  function 
of  time,  related  to  ^  e(r,t)  by  a  transformation  of  the  type  (U).  Explicitly  we 

have 

I(t)  =  l'(t)   +  i  I*  (t), 

where  l'(t)  and  I  (t)  are  the  real  and  imaginary  parts  of  l(t).  From  the  fact 
that  <  A  e(r,t)  >  ■=  0,  it  follows  that  <  l(t)  >  »  0.  On  the  other  hand,  the 
average  envelope  is 


<  |l(t)|  >  '  </l'^(t)  +  l"^(t)  >  f     0. 

In  the  opinion  of  the  author,  none  of  the  many  derivations  of  equation  (16) 
found  in  the  radio  scattering  literature  is  sufficiently  free  of  heuristics  to 
be  considered  definitive.  However,  for  the  purposes  of  this  paper,  we  shall  re- 
gard the  familiar  formula  (l6)  as  a  correct  result.  Two  of  the  many  assumptions 
on  which  the  validity  of  the  approximation  (16)  rests  are  worth  noting.  One  is 
that  the  maximum  scattering  angle  is  small,  which  permits  the  neglect  of  polar- 
ization effects.  The  other  is  that  the  spatial  pattern  of  ^  e(r,t)  remains  sub- 
stantially unchanged  in  the  longest  time  it  takes  the  incident  radiation  to 
traverse  the  scattering  volume  V,  which  permits  the  neglect  of  retardation  effects. 
We  now  describe  two  further  approximations  that  can  be  applied  to  (l6). 

For  an  arbitrary  (sure)  function  h(r),  consider  the  quantity 

5h  =  Max   h(r)  -  h(r')  ,  r,  r'  in  V, 

Then  if  6h  «  h(r  ),  where  r  is  the  center  of  V,  it  is  a  good  approximation  to 
c        c 

replace  h(r)  by  the  constant  value  h(r  )  in  any  integral  over  V,  In  particular, 
if  the  inequalities  6R  «  R(r^),  6r'  «  ^'(rp),  and  6K  «  K(r^)  are  valid,  the 
functions  R(r),  R'(r)  and  K(r)  can  be  replaced  by  their  values  at  r^  in  equation 
(16),  an  approximation  which  ve  shall  call  the  Fraunhofer  approximation,  for  an 
evident  reason.  The  Fraunhofer  approximation  is  a  good  approximation  in  some 
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radio  scattering  experiments  and  not  in  others,  as  the  following  simple  consider- 
ations show. 

Consider  a  scattering  experiment  wi.th  path  length  D  and  with  identical  antennas 
of  half -power  beamwidth  p.  Then  R(r  )  -^  R  (r  )  ~  D/2,  and  5R  ~  5R  ~  (np/2  @  .  ), 

so  that  the  conditions  5R  «  R(r  )  and  6R  «  R  (r  )  become  S  «  (S)  .  ,  Moreover. 

c  c  '^  min  ' 

(H)(r  )  -  O    .      +  B,   and  6  (£)  '-  2p,    so  that  the   conditions  5  (S>  «  ®(r  )   also  re- 
c  mm       ^ '  '  '  ^   c ' 

duces  to  the  inequality  B  «  <3)    .    ,     As  ve  have  noted,  (H)     •      ~  D/a,   so  that  this 

"^      rain  '    mm 

inequality  is  reasonably  valid  for  long  distance  experiments,  using  high  gain 

antennas.  However,  for  short  distance  tropospheric  experiments,  the  inequality 

may  fail  completely,  e.g.,  if  p  ^   l/5o  radian  and  D  •^    100  km. 

Another  common  approximation,  which  we  shall  call  the  uniform  illumination 

approximation,  for  an  evident  reason,  consists  in  settinp  F(r)  «■  1,  throughout 
V.  The  Fraunhofer  and  uniform  illumination  approximations  can  never  be  equally 
good.  If  the  Fraunhofer  approximation  is  good,  there  is  considerable  variation 
of  F(r)  throughout  V,  by  the  very  definition  of  V.  On  the  other  hand,  if  the 
Fraunhofer  approximation  fails,  the  variation  of  lf(r)  id.thin  V  determines  a  sub- 
volume  of  V  which  contributes  most  of  the  scattered  radiation,  and  then  the 
approximation  F(r)  '^  1  is  pood  in  this  effective  subvolume  ^  -' , 


6.   Statistical  structure  of  the  scattered  signal 

We  turn  now  to  a  study  of  the  statistical  structure  of  l(t).  Of  basic  im- 
portance is  the  fact  that  l(t)  is  a  gaussian  random  function  of  timej  this  is 
shown  as  follows » 

Consider  a  radio  scattering  experiment  with  path  length  D,  identical  antennas 
of  half -power  beamwidth  p,  and  typical  scattering  angle  C©  .  Then  the  scattering 
volume  V  can  be  estimated  as 
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Let  N  be  the  ratio  of  V  to  a  correlation  volume,  i.e., 

N  -  VA^  -  (l/8®)(DpAo)^  • 

In  any  ratio  scattering  experiment  N  »  1;  for  example,  if  p  —'  lAOO  radian  and 
D  ~  200  km,  then  N  ~  10  .  We  now  divide  V  into  subvolumes  V^,.<..,V  j  for  simpli- 
city and  vdthout  loss  of  generality,  we  imagine  all  these  subvolumes  to  be  cubes 

of  side  L  .  Let  the  contribution  to  l(t)  from  the  volume  V  be  denoted  by  I. (t), 
o  ^  •'■ 

so  that 

N 

I(t)  -  zi     h^^^    • 

I.(t)  is,  of  course,  given  by  equation  (16)  with  the  volume  of  integration  replaced 

by  v.,  i'or  every  fixed  t,  the  random  variable  I.(t)  is  substantially  independent 

of  the  random  variable  I.(t),  j  /  i,  because  of  our  choice  of  the  volumes  V^,  uraess 

V  and  V  are  adiacent,  in  which  case  there  is  a  weak  dependence  between  I. (t)  and 
i     0  ^ 

I.(t).  It  follows  by  the  central  limit  argument  that,  for  every  t,  l(t)  is  approxi- 
mately  gaussian.  Indeed,  by  the  n-dimensional  central  limit  argument,  the  joint 
distribution  of  the  random  variables  l(t^),.. .,l(t^) ,  where  n  «  N  and  t^,...,t^ 
are  any  n  instants  of  time,  is  approximately  gaussian,  i.e.,  l(t)  is  approximately 
a  gaussian  random  function  of  time  .  The  received  scattered  signal  looks  like 
narrow-band  gaussian  noise,  and  numerous  experiments  confirm  the  fact  that  the 
univariate  distribution  of  l(t)  is  gaussian,  i.e.,  that  the  univariate  distribution 
of  ll(t)l  is  Rayleigh"-  -'*'-^^-'}  the  author  knows  of  no  experiment  that  tests  the 
multivariate  gaussian  nature  of  l(t). 


*By  the  central  limit  argument  we  mean  the  qualitative  statement  that  the  sum  of 
a  large  nioiriber  N  of  random  variables  is  approximately  gaussian  if  the  dependence 
between  the  summands  is  sufficiently  weak  and  if  no  small  group  of  summands  con- 
tributes most  of  the  variance  of  the  sum.  (in  the  application  given  here,  we 
assume  that  the  weighting  factors  F(r),  R(r),  R  (r) ,  etc.,  do  not  vary  enough  to 
inhibit  the  operation  of  the  central  limit  argument.)  This  rough  statement  is 
supported  by  a  highly  developed  central  limit  theory,  the  sophistication  of  which 
far  transcends  the  ordinary  needs  of  physicists  and  engineers.  By  the  n-dimensional 
central  limit  argioment,  we  mean  the  extension  of  the  ordinary  central  limit  argument 
to  the  case  of  n-dimensional  A'ectors,  where  n  «  N» 
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The  eddy  heuristics  that  lead  us  to  believe  that  dielectric  turbulence  has 
the  local  harmonic  property  suggest  that  the  subvolumes  used  in  the  central  limit 
argument  need  only  have  linear  dimension  d,  where  K   «  d  «  L  .  (In  most 
scattering  experiments,  such  values  of  d  can  be  found.)  If  so,  this  would  result 
in  more  summands  and  a  more  povrerful  central  limit  argument.  However,  because 
of  the  usual  large  size  of  N,  the  cruder  central  limit  argument,  which  uses 
correlation  volumes,  is  quite  satisfactory. 

In  some  earlier  workL  ■'*'-  -',  the  author  presented  an  argument  purporting 
to  show  that  l(t)  is  approximately  gaussian  if  the  dielectric  noise  originates 
in  turbulent  pressure  fluctuations;  the  central  limit  argument  was  used  to  show 
that  Fourier  coefficients  of  the  turbulent  pressure  field  are  gaussian  random 
variables.  Hovrever,  the  validity  of  this  argument  seems  to  be  open  to  question  . 
In  any  event,  this  matter,  although  of  some  interest  in  turbulence  theory,  is 
of  only  academic  interest  in  radio  scattering  theory,  since  there  is  now  little 
doubt  that  the  turbulent  pressure  field  makes  a  negligible  contribution  to  atmos- 
pheric dielectric  noise. 

Since  l(t)  is  a  zero-mean  gaussian  random  function,  then,  as  is  well  known, 
its  statistical  structure  is  completely  specified  by  giving  its  covarianca 
<  l(t)I*(t  )  >.  Moreover,  since  l(t)  is  a  linear  transform  of  ^  e(r,t),  we 
need  only  know  the  covariance  (or  spectrum)  of  ^  £(r,t)  in  order  to  calculate 

Fourier  coefficients  of  the  pressure  field  can  be  written  as  convolutory  sums  of 
the  Foiirier  coefficients  of  the  turbulent  velocity  field  in  a  familiar  manner*-  -^'^    -^. 
It  is  correctly  observed  in  (Y]  that  the  velocity  coefficients  are  not  statisti- 
cally independent  as  assumed  in  [2?] .  However,  it  is  very  difficult  to  ascertain 
whether  the  dependence  between  the  summands  is  weak  enough  to  meet  the  require- 
ments of  the  dependent  central  limit  argument,  as  assumed  in  [5] .  The  fact  that 
the  univariate  distribution  of  the  velocity  field  is  gaussian,  whereas  its  bi- 
variate  distribution  at  small  separations  is  non- gaussian,  shows  how  delicately 
poised  such  arguments  are. 
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the  covariance  (or  spectrum)  of  l(t).  Two  consequences  of  this  fact  have  already 
been  mentioned,  nainely  that  all  station?rity  and  ergodicity  properties  of  ^  e(r,t) 
need  only  be  second-order,  and,  secondly,  that  A  e(r,t)  can  be  replaced  by  any 
process  with  the  same  covariance  without  affecting  the  scattered  signal;  in  parti- 
cular, we  can  substitute  the  zerc-mean  gaussian  process  with  the  same  covariance, 
which  always  exiFts  (see  [6],  p.  7U)» 

To  complete  our  study  of  the  statistical  structure  of  the  scattered  signal, 
we  new  relate  the  covariance  <  l(t)l  (t  )  >  to  the  ndxed  spectral  density  of  the 
process  A  e(r,t),  defined  in  section  U«  We  proceed  in  two  steps.  First  we  make 
the  Fraunhofer  and  uniform  illumination  approximations  by  setting  Fir)  =  1, 
R(r)  =  R,  R  (r)  -  R  ,  and  K(r^  =  ff,  where  by  R,  R  ,  and  K,  we  mean  the  values 
of  the  corresponding  functions  at  r  ,  the  center  of  the  scattering  volume.   (We 
assume  that  the  scattering  volinne  is  a  cube  of  side  d,  centered  at  r  ;  clearly 
this  entails  no  loss  of  generality.)  We  also  assume  that  ^  e(r,t)  is  locally 
stationary  in  the  scattering  volume  V,*  so  that  it  can  be  replaced  in  V  by  a 
"  piece"  of  a  globally  stationary  random  function  with  the  same  covariance. 
(From  now  on  it  is  understood  that  all  stationarity  requirements  need  only  apply 
in  the  second-order  sense.)  ^'ith  these  assumptions  it  follows  from  (l5)  and  (16) 
that 

(17)  <  I(t)l*(t')  >  '-  (tiDA^Rr')^  <  J(d,K;t)J*(d,Kjt')  > 

~ (nD/X^RR' )^   (2n)^  V  <  ^  e^  >  |(K, t ), 

where  X  =   t-t  . 

Next  we  consider  the  case  where  the  functions  F(r),  R(r),  R  (r),  K(r)  and 
the  dielectric  noise  vary  appreciably  in  the  scattering  volume  V.  For  simp.Licity, 
vre  confine  ourselves  to  a  study  of  the  average  power  <|l|  >.  As  before  we 
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divide  V  into  N  subcubes  of  side  L^.  It  is  easily  seen  that  in  volumes  of  side 
L^-~'3CXD  meters,  there  is  negligible  percentage  variation  of  the  functions  F(t') , 
R(r),  R  (r),  and  K(r),  regardless  of  whether  the  Fraunhofer  or  uniform  illumination 
approximation  applies  in  the  total  voltune  V.  Thus,  in  each  of  the  subvolumes  V. , 
the  functions  F(r),  R(r),  R  (r),  and  K(r)  can  be  replaced  by  their  values  at  the 
center  of  V^,  which  we  denote  by  F^,  R^,  R^,  and  K. .  It  follows  that 

N 
(18)        <|T|2>  ^  (nD/X^)^  ^  (^i/\\)^  "   Ji(Lo,K^^t)J^(L^,K^jt)>  , 

i=l 

where  the  subscript  on  J.    reminds  us  that  it  pertains  to  the  dielectric  noise  in 
the  subvolume  V..    (We  neglect  interactions  between  different  J..) 

To  reduce   (l8)  further,  we  recall  that  I^  is  the  least   linear  dimension  of  the 
domains  of  local  stationarity  of  ^  e(r,t)   as  well  as  its  correlation  distance,   so 
that  the  volumes  V     are  certainly  domains  of  local  stationarity  of  A  e(r,t).     (Our 
subdivision  of  V  is  deliberately  conservative!  in  most  scattering  experiments,   the 
domains  of  local  stationarity  of  ^  c(r,t)   can  be  expected  to  be  ruch  larger.)     Thus 
/\  £(r,t)  can  be   approximated  in  each  of  the  volumes  V.   by  a  piece  of  a  suitable 
globally  stationary  process.     We  now  assume  that  for  our  purposes  the  change  in 
dielectric  noise  from  subvolume   to  subvolume  can  be  regarded  as  due  only  to  a  change 
of  the  noise  power  <  ^  e   (r)  >,  and  not  to  a  change  in  the  normalized  spectral 
density  J  (k,0).       Then,    noting  that  in  most  radio  scattering   experiments 
KL^/2n  ^  €>L^/X  »  1,   where     (3)    is  a  typical  scattering  angle,   and  assuming  that 
dielectric  turbulence  has  the  local  haimonic  property  (see  Sect.   U),    we  can  apply 
the  arriroxiraation  (1$)  to  the  covariance  of  J.,    so  that   (l8)  becOTies 


(19) 


N 


<|ir  >  r-   (nDA^)^(2n)^  ^I  V.  (F^/R^R.)^  <  A  ^  >i  I  (^.,0), 


2  ?  ^ 

where  <  /^t     >.  denotes  the  value  of  <  ^e  (r)  >  at  the  center  of  V.,  Equation 

(19)  can  be  approximated  by  the  integral 

The   plausibility  of  this  assumption  rests  in  part  on  the  fact  that  we  shall  only 
require  the  spectrum  of  ^^  E(r,t)  to  be  defined  at  wave  numbers  K  satisfying 

KL  »  2n. 

0 
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(20)  <|I|2>     ~  (.DA')^(2n)^  ^'^!  ^    <  /^,^{t)>  l[K{r),0]   dr. 

{    R(i^R   (r)        "^  ^L  J 

If  we  define  an  effective   scattering  volume  V  _„  by 


rr' 


«^f       i^  R(?)R   (i^       <A^>  i(K,0) 


where  <  /\  e     >  is  now  the  value  of<^e(r)>atr,    then,   >rtien    2^  "  0,    (20)  Is 
just   (17)  with  V  replaced  by  V  „_.     The  usual  rovigh  estimates  of  the  scattered 
power  correspond  to  the  approximation  V  „_  "^  V. 


7 .       Conclusions 

From  the  mathematical  standpoint,  much  of  radio  scattering  theory  consists 
in  a  study  of  finite,  vreighted  Fourier  space-transforms  of  a  certain  random 
function  of  space   and  time,  the  dielectric  noise.     This  transform  is,  of  course, 
itself  a  random  function  of  time,  the  scattered  signal.     Our  purpose  has  been  to 
f?ive   a  unified  treatment  of  the   relation  between  these  two  random  functions,  ex- 
pressed in  the  appropriate  mathematical  language.     At  points  we   have  had  to  sub- 
stitute plausible  heuristics  for  rigor,  becaus-  no  theory  of  the  local  properties 
of  second-order  random  functions,  in  particular  local  stationarity,  and  what  we 
call  the  local  harmonic  property,  has  yet  been  developed.     The  complicated  radio 
scattering  problem  is  obviously  not  the  context  in  which  to  attempt  these  exten- 
sions of  random  function  theory. 
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ATTN:    Tech.   Library,  Delmer  C.    Ports 

Henry  Jasik,  Consulting  Engineerini; 
61  Second  Street 
i'dneola.  New  York 

Jet  Propulsion  Laboratory 
I48OO  Oak  Grove  Drive 
Pasadena,  California 
ATTN:     I.E.  Newlan 

Martin  Katzin  and  Co.,   Consultants 
711     ILth  Street,  N.W. 
Washington  5,  D.C. 

Lincoln  Laboratory 

P.O.   Box  73 

Lexington  73,  Massachusetts 

ATTO:   Henry  Straus,  Document  Rm.,  A-229 

Lockheed  Aircraft  Corporation 
.'•Bsaile  Systems  Division 
7701  Woodley  Avenue 
Van  Nuys,  California 
ATTN:   M3D  Library  -902,  Dept.   82-10 
E.A.  Blasi,  Head,  Antenna  Lab. 

Lockheed  Aircraft  Corporation 

Call  fomia  Division 

Engineering  Library 

Dept.  72-25,  Plant  A-1,  Bldg.  63-1 

ATTN:  N.C.  Harnoia 

The  Glenn  L.  Martin  Company 

Denver,  Colorado 

ATTN:  Wm.  V.  Foley,  Engineering  lib. 

The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
ATTN:  Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Mfg.  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
ATTNt  H.  Warren  Cooper 

Mathematical  Reviews 
80  Waterman  Street 
Providence  6,  Rhode  Island 

The  W.L.  Maxon  Corporation 
160  West  31ith  Street 
New  York,  New  York 
ATTN:  Antenna  Laboratory 

McDonnell  Aircraft  Corporation 

Box  516 

St.  Louis  3,  Missouri 

ATTN:     R.D.  Detrich 

Engineering  Library 

McMillan  Laboratory,   Inc. 
Ipswich,  Massachusetts 
"ATTN:     Security  Officer 
Docume'^t  Room 

Melpar,  Inc. 

3000  Arlington  Boulevard 

l^lls  Church,   Virginia 

ATTN:  Engineering,  Technical  Lib. 


Micrbwave  Radiation  Company 
19223  So.   Hamilton  Street 
Gardena,   California 
ATTN;     Mr.   Morris   J.  Ehrlich,   Pres. 

Northrop  Aircraft,   Inc. 
Engineering^  Department 
1000  E.  Broadway 
Hawthorne,  California 
ATTN:     E.A.  Preitas,  Engineering 
Library,   Plant  III 

North  American  Aviation,  Inc. 
12?li4  Lakewood  Boulevard 
Downey,  California 
ATTtl:      Technical   Library, 
Dept.   991-h 

North  American  Aviation,   Inc. 
Los  Angeles  International  Airport 
Los  Angeles  U5,  California 
ATTN:     Engineering  Technical  File 

Page  Communications  Engineers 
71c  lljth  Street,   N.W. 
Washinrton  5,  D.C. 
ATTN :      Library 

Philco  Corporation  Research  Division 
A  and  Lippincott  Street 
Philadelphia  31),  Pa. 
ATTO:     Technical  Library 

Pickard  and  Bums,   Inc. 

21jO  Highland  Avenue 

Needhaji'.  91),  Mass. 

AT  El:     Dr.   J.T.   de  Bettencourt 

Polytechnic  Research  and  Development 

Company ,   Inc. 

202   Tiilary  Street 

Brooklj'n  1,  New  York 

ATTO:     Technical  Library 

Radiation  Engineering  Laboratory 
Maynard,  I'lass. 
ATTN;     Dr.   John  Ruze 

Radiation,    Inc. 
P.O.   Drawer  37 
Melbourne,  Florida 
ATTN:      Technical  Library 

Mr.   M.L,   Cox,  Assistant 

Project  Engineer 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  Ycrk 

ATTN:     P.S.   Carter,  lab.   Library 

RCA  Laboratories 
Princeton,   New  Jersey 
ATTN:     Miss  Fern  Cloak,   Librarian 
Research  Library 

RCA  Defense  Electronic  Products 
Camden  2,  New  Jersey 
ATTN:     Miss   J.K.   Steever,   Librarian 
Library  lC-2 

(2)rhe  Rair\o-Woolfirid'^e  Corporation 
S730  Arlior  Vitae  Street 
Los  Angeles  kSj   California 
ATTN:     Margaret  C.  Whitnah, 

Chief  Librarian,  R^^search 
and  Development  Library 

The  Rand  Corporation 

Via:     SBAMA  Liaison  Office 

1700  Main  Street 

Santa  Monica,  California 

ATTU:     Library 

Rantic  Corporation 

Carlatasos,   California 
ATTN:      Technical  Library 

Raytheon  Manufacturir.g  Company 
M  and  R  Division 
Bedford,  Massachusetts 
ATTN:      Irving  Goldstein 

Raytheon  Laboratory 
Wayland,  Massachusetts 
ATOI:     Document  Lib.,    John  E.  Walsh 
Dept.   8971),  Antenna 
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Raytheon  Manufa.-turing  Company 
Wayland  Laboratory 
Wayland,   Massachusetts 
ATTN:     Miss  Alice  G.  Anderson,   Lib. 
Wayland  Library 

Republic  Aviation  Corjioration 

FarmingdalG,   L.I.,  New  York 

ATTN:     H.   Stieglitz,  Engineering  Lib. 

Rheem  Manufacturin^r  Company 
Researc:;  and  Development 
9236  East  Hall  Road 
Downey,   California 
ATTN:     J.C.   Joerger 

Ryan  Aeronautical  Company 
San  Diego,   California 
ATTN:      Technical  Library 

Sanders  Associates 
137  Canal  Street 
Nashua,   New  Hampshire 
ATTN:     N.R.  Wild,   Library 

Sandia  Coritoration 

P.O.  Box  5800 

Classified  Doojment  Division 

Albuquerque,  New  Mexico 

Sperry  Gyroscope  Company 
Great  Neck,  L.I.,  New  York 
ATTN:  Engineering  Library 

Florence  W.  Turnbull, 

Engr.  Librarian 

StarLford  Research  Institute 
Menlo  Park,  California 
AT1TJ:  Engineering  Document  Library 
Mrs.  Beatrice  Gibson,  Libr. 

Sylvania  Electric  Company 
100  First  Street 
Walt ham,  Mass. 
ATT:  Librarian 

Systems  Laboratories  Corporation 
150I6  Ventura  Blvd. 
Sherman  Oaks,  California 
ATTN:  Dr.  W.C.  Hoffman 

Technical  Research  uroup 
17  Union  Square  West 
New  York  3,  New  York 

A.S.  Thomas,  Inc. 
161  Devonshire  Street 
Boston,  Mass. 
ATTN:  A.S.  Thomas 

Westinghouse  Electric  Corp. 
2519  Wilkens  Avenue 
Baltimore  3,  Maryland 
ATTN:  Engineering  Library 
Kent  M.  Mack 

Wheeler  laboratories.  Inc. 
122  Cutter  Mill  Road 
Great  Neck,  New  York 
AWN:  Mr.  Harold  A.  Wheeler 

Zenith  Plastics  Company 
Box  91 

Gardena,  California 
ATTN:  Mr.  S.S.  Oleesky 

Library 

Geophysical  Institute 
University  of  Alaska 
College,  Alaska 

Electronic  Research  Laboratory 
Div.  of  Electronic  Engineering 
University  of  California 
Berkeley  h,   California 
ATTN:  S.  Silver 

Electronics  Research  Laboratory 

232  Cory  Hall 

University  of  California 

Berkeley  U,   California 

ATTN:  J.R.  Whinnerv 

California  Institute  of  Technology 

1201  E.  California  Street 

Pasadena,  California 

ATTN:  Dr.  Charles  H.  Papas 


Carnegie  Institute  of  Technology 
PittsViurgh  13,   Pennsylvania 
ATTN:     Albert  E.   Heins 

School  of  Electrical  Engineering 

Cornell  University 

Ithaca,   New  York 

AT  IN:     Dr.   Henry  G.  Booker 

College  of  Engineering 
University  of  Florida 
Gainesville,   Florida 
ATTN:     Engineering  Sciences  Library 
Professor  M.H.   Latour 

Engineering  Experiment  Station 
jeorgia  Institute  of  Technology 
Atlanta,  Georgia 
ATTN:       Georgia  Tech.   Library 
Mrs.   J.H.  Crosland 

Technical  Reports  Collection 

Har'/ard  University 

Gordon  l>1cKay  Library 

303A  Pierce  Hall 

Oxford  Street 

Cambridge  38,   Massachusetts 

ATTO:      Mrs.  E.L.   Hufschmidt,   Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge,   Massachusetts 
ATTN:     Dr.   Fred  L.  Whipple 

Serials  Department     -220-S  Library 

University  of  Illinois 

Urbana,   Illinois 

ATTN:     Engineering  Library 

Antenna  Section 

Electrical  Engr.  Research  laboratory 
University  of  Illinois 
Urbana,   Illinois 
ATTO:     EERL  Ubrary,  Room  210, 
Dr.  R.H.  DuHamel 

The  Johns  Hopkins  University 
Department  of  Physics 
Baltimore  18,   Maryland 
ATTN:     Professor  Donald  E.  Kerr 

Radiation  Laboratory 
Johns  Hopkins  University 
1315  St.   Paul  Street 
Baltimore  2,  Maryland 

Applied  Physics  laboratory 
The  Johns   Hopkins  University 
8621  Georgia  Aven'.e 
Silver  Spring,  i-iaryland 
ATTN:     Document  Library 

Massachusetts  Institute  of  Technology 
Research  Laboratory  of  Electronics 
Document  Room  20B-221 
Cambridge  39,  Massachusetts 

University  of  I'tlchigan 
Electronic  Defense   Group 
Engineerinj^  Research  Institute 
Ann  Arbor,  Michigan 
ATTN;     J. A.  Boyd,   Supervisor 

Willow  Run  labs.,   of  the  Engineering 
Research  Institute 
The  University  of  Michigan 
Willow  Run  Airport 
ipsilanti,  Miichigan 
ATTN:     K.M.   iie.cel.   Head,   Theory 
and  Analysis  Department 

University  of  Michigan 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilanti,   Michigan 
ATTN:     Shelia  Coon,   Librarian 

University  of  Minnesota 
Minneapolis,  Minnesota 
Engineering  Library 
ATTN:      Head   Librarian 
Microwave  laboratory 
Electrical  Engineering  Dept. 
Northwestern  University 
Evanston,   Illinois 
ATTN:     Prof.  R.E.  Beam 
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Antenna  Laboratory 

Department  of  Electrical  Engineering 

Ohio  State  University 

Colujnbus,  Ohio 

ATTN:     Dr.  T.E.  Tice 

Oklahoma  University  Research  Found. 

Norman,  Oklahoma 

ATTNi  Technical  Library 

Professor  C.L.  Farrar 

Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,  New  Yrrk 
ilTN:  Dr..  A. A.  Oliner 

Polytechnic  Institute  of  .Brooklyn 

55  Johnson  Street 

Brooklyn  1,  New  York 

ATIN:  Microwave  Research  Institute 

Syracuse  University 
Dept.  of  Klectrical  Engr. 
Syracuse  10,  New  York 
ATTN:  Technical  Library 

University  of  Texas 
Electrical  Engr.  Research  Lab. 
P.O.  Box  8026,  Univ.   Station 
Austin  12,    Ttexas 
ATTO:     John  R.   Gerhardt 

Dept.  of  Physics 
University  of  Texas 
Austin  12,  Texas 
ATTN:     Physics  Library 
Claude  W.  Morton 

Tufts  University 

Hedford  55,  Massachusetts 

ATTN:     Research  Lab.   of  Physical 

Electronics 

Professor  Charles  R.   Hingins 

University  of  Toronto 
Dept.   of  Electrical  Engineerinp, 
Toronto,   Ontario,   Canada 
ATTN:   Prof.   George  Sinclair 

University  of  Washington 
Dept.  of  Electrical  Engineering 
Seattle,  Washington 
ATM:      Prof.   0.  Held 

Ionosphere  Research  Laboratory 
Pennsylvania  State  Collage 
State  College,   Pennsylvania 

ATTN:   Prof.  A.H.  Waynlck,  Director 

Institute  of  Mathei>iatical   Science 

25  Waverly  Place 

New  York  3,  New  York 

ATTN:     Ib-a.   JoAn  Segal,   Librarian 

School  of  Electrical  Engineering 

Purdue  University 

iafayette,  Indiana 

ATTN:  Professor  F.V.  Schultz 

Department  of  Electrical  Engineering 
California  Institute  of  Teclinology 
ATTN:  Dr.  Charles  H.  Papas 

Electronics  Division 
Rand  Corf^oration 
1700  Main  Street 
Santa  Monica,  California 
ATTN:  Dr.  Robert  Kalaba 

National  Bureau  of  Stai  dards 

Washington,  D.C. 

ATTN:  Dr.  W.K.  Saunders 

Applied  Mathematics  and  Statistics  Lab. 

Stanford  University 

Stanford,  California 

ATTN:  Dr.  Albert  H.  Bowker 

Department  of  Physics  and  Astronomj' 
Michigan  State  University 
East  lansinr;,  Michigan 
ATTO:     Dr.  A.  Leitner 


University  of  Tennessee 
Knoxville,  Tennessee 
ATit):     Dr.  Fred  A.  Tlcken 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
ATTN:     Dr.  A.  Erdelyl 

Wayne  University 

Detroit,  Michigan 

ATTN:     Professor  A.K.   Stevenson 

Mathematics  Department 
Stanford  University 
Stanford,   California 
ATTO:     Dr.    Harold  Levine 

University  of  Minnesota 
Minneapolis  llj,   Minnesota 
ATTN:   Prof.   Paul  C.  Rosenbloora 

Department  of  Mathematics 
University  of  Pennsylvania 
Philadelphia  h.   Pa. 
ATTN:     Professor  Bernard  Epstein 

Applied  Physics  laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,   Maryland 
ATTN:     Dr.  B.S.  Gouray 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,  D.C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Tech. 
Cambridge   39,   t^ssachusetts 
ATTN:     Dr.   L.J.  Chu 

Nuclear  Development  Assoc.   Inc. 

5  New  Street 

White  Plains,  New  York 

ATTN:   Library 

Lebanon  Valley  College 
Annville,  Pennsylvania 
ATTN:  Professor  B.H.  BissiniTer 

Dept.  of  Physics,  Tliaw  Hall 
University  of  Pittsburgh 
Pittsburgh  13,  Pa. 
ATTN:  Dr.  Edward  Gerjuoy 

Dept.  of  Physics 
Amherst  College 
Amherst,  Massachusetts 
ATTN:  Dr.  Arnold  Arons 

California  Institute  of  Tech. 
Electrical  Engineering 
Pasadena,  California 
AT  IN:  Dr.  Zohrab  A.  Kaprielian 

Dr.  Rodman  Doll 
209A  Emmet  St. 
Tpsilanti,  Michigan 

California  Institute  of  Tech. 
Pasadena  U,   California 
ATTN:  Mr.  Calvin  Wilcox 

(3) Mr.  Robrrt  Brockhurst 

Woods  Hole  Cceanographic  Inst. 
Woods  Hole,  Massachusetts 

National  Bureau  of  -"tandards 
Boulder,  Colorado 
ATTN:  Dr.  R.  Gallet 

Engineering  Library 
Plant  5 

Grumman  Aircraft  Corp. 
Bethpage,  L.I.,  New  York 
ATTN:  Mrs.  A.M.  Gray 

Mrs.  Jane  Scanlon 
281:  South  Street 
Southbridge,  Massachusetts 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  California 


Microwave  Laboratory 
55  Johnson  Street 
Brooklyn,   New  York 

ATTN:     Dr.  Bernard  Llppraann 

University  of  Minnesota 
The  University  Library 
Minneapolis  11:,  Minnesota 
ATTO:     Exchang3  Division 

Professor  Bernard  Friedman 
55  Hilltop  Avenue 
New  Roc hell e,   NewYork 

Lincoln  Laboratory 

liassachusetts  Institute  of  Tech. 

P.O.   Box  73 

Lexini^on  73,  Mass. 

A'JTN:     Dr.  Shou  Chin  Wang,     Rm.  C-351 

Melpar,  Inc. 

3000  Arlington  Boulevard 

Falls  Church,    Virginia 

ATTN:     Mr.  K.S.   Kelleher,   Section  Head 

Antenna  laboratory,  ERD 

HQ  Air  Force  Cambridge  Research  Center 

Laurence  C.   lianscom  Field 

Bedford,  Massachusetts 

ATTN:     Nelson  A.   Logan 

Electronics  Research  Directorate 

HO  Air  Force  Cambridge  Research  Center 

I^uTence  0.   Hanscon  Field 

Bedford,  Massachusetts 

ATTN:     Dr.   Philip  Newman,   CRRK 

HQ  Air  Force  Cajnbridge  Research  Center 
Laurence  C.   iianscom  Field 
Bedford,  f-Iassachusetts 
ATTN:     Mr.   fVancis  J.   Zucker 

Crosley  AVCO  Research  labs. 
257  Orescent  Street 
Waltham,   Massachusetts 
ATTN:     Mr.  N.C.   I'lerson 

Antenna  Research  Section 
Microwave  Laboratory 
Hughes  Aircraft  Company 
Culver  City,  California 
ATTO:     Dr.  Richard  B.  Barrar 

Columbia  University 
Hudson  Laboratories 
P.O.  Box  239 
lliS  Palisade  Street 
Dobbs  Kerry,  New  York 
ATTN:     Dr.   M.W.   Johnson 

Institute  of  Fluid  Dynamics  and  Applied  Hath. 
University  of  Maryland 
College  Park,  Maryland 
ATTN:     Dr.  Elliott  tontroU 

Dept.   of  Electrical  Engineerini; 
Wasliington  University 
Saint  Louis  5,   Missouri 
ATTN:     Professor  J.   Van  Bladel 

Dept.   of  the  Navy 

Office  of  Maval  Research  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,  California 

Brandeis  University 
Waltham,   Massachusetts 
at™  :     Library 

General  Electric  Company 
Microwave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,   California 
AT  IN:      Library 

Hughes  Research  Laboratories 
Hughes  Aircraft  Conpany 
Culver  City,   California 
ATTN:     Dr.  W.A.  Dolid 
Bldg.   12,     Rm.  2529 

Sr^yth  Research  Associates 
3930  lith  Avenue 
.San  Diego  3,   California 
ATTN:     Dr.   John  B.  Smyth 


Electrical  Eneineeriru; 
California   Institute  oi    ici^i 
Pasadena,   California 
ATTN:     Dr.  Georges  G.  Weill 

Naval  Research  Laboratory 
Washini^ton  25,  D.C. 
ATTN:     Dr.  Henry  J.   Passerini 
Code  5278  A 

Dr.  George  Kear 
10585  N.   Stelling  Rd. 
Cupertino,   California 

Brooklyn-  Polytechnic   Institute 
85  Livingston  Street 
Brooklyn,  New  lork 
ATTN:     Dr.  Matlian  Marcuvitz 

Department  of  Electrical  Engp. 
Brooklyn  Polytechnic   Institute 
85  Livingstcn  Street 
Brooklyn,    New  York 
AT  TO:     Dr.    Jerry  Shraoys 

Dept.  of  Itetheraatics 
University  of  New  >fexico 
Albuquerque,  New  Mexico 
ArTN;     Dr.   I.  Kolodner 

W.L.  ffaxon 

I46O  W.  3lith  Street 

New  York,  N.Y, 

ATTO:     Dr.    1-krry  Hochstadt 
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